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ABSTRACT 

Hyper- accretion discs around black holes emit copious neutrinos and anti-neutrinos. 
A fraction of the emitted neutrinos convert to electron-positron plasma above the disc 
through the annihilation reaction vV e~^e~ . This process may drive relativistic jets 
associated with GRB explosions. We calculate the efficiency of energy deposition by 
neutrinos. Our calculation is fully relativistic and based on a geodesic-tracing method. 
We find that the efficiency of neutrino heating is a well-defined function of (i) accretion 
rate and (ii) spin of the black hole. It is practically independent of the details of 
neutrino transport in the opaque zone of the disc. The results help identify accretion 
discs whose neutrino emission can power GRBs. 

Key words: accretion discs, black hole physics, gamma-ray burst: general magnetic 
fields, neutrinos, relativistic processes 



1 INTRODUCTION 

A plausible model for the central engines of gamma-ray 
bursts (GRBs) pictures a transient, hyper- accreting disc 
formed around a rotating black hole (see Beloborodov 2008 
for a recent review) . The disc is a source of copious neutrinos 
and anti-neutrinos, which partially annihilate above the disc 
and turn into pairs, u -\- u ^ e~ -\- . This process was 
proposed as a possible mechanism for creating relativistic, 
-dominated jets that could power observed GRBs (Eich- 
ler et al. 1989). However, its efficiency remained unsettled. 
Its accurate calculation requires a detailed relativistic model 
for the neutrino source — the accretion disc — as well as 
tracing the neutrino propagation in the Kerr spacetime of 
the black hole. 

A detailed relativistic model for GRB discs was com- 
pleted recently (Chen & Beloborodov 2007, hereafter CB07). 
It describes the disc down to its inner edge and gives accu- 
rate energy fluxes carried away by v and u at all radii r. In 
the present work we trace the neutrino trajectories and cal- 
culate the rate of uiy annihilation around the disc. Neutrino 
annihilation was previously calculated in a number of works 
(e.g. Popham, Woosley & Fryer 1999.; Asano & Fukuyama 
2001; Birkl et. al. 2007). Our work has three motivations: (i) 
A relativistic calculation of uiy annihilation has never been 
done for a realistic accretion disc around a spinning black 
hole. Previous works either used a toy model for neutrino 
source (e.g. Birkl et al. 2007) or replaced neutrino trajec- 
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tories by straight lines (Popham et. al. 1999). (ii) The effi- 
ciency of lyjy annihilation strongly depends on the accretion 
rate M and the spin parameter a of the black hole. It is de- 
sirable to know this dependence and identify the range of M 
and a where lyjy annihilation can provide the observed en- 
ergy of GRB explosions, (iii) In addition to uiy annihilation, 
neutrinos can create pairs off magnetic field (Kuznetsov 
& Mikheev 1997; Gvozdev & Ognev 2001). The contribu- 
tion of this process to the energy deposition rate should be 
included. 

The paper is organized as follows. Section |2] describes 
the setup and method of our calculations. The results are 
presented in Section [3] and summarized in Section 4. 



2 MODEL DESCRIPTION 

Neutrinos emitted by the disc follow geodesies in Kerr space- 
time. The efficiency e of their annihilation can be calculated 
numerically by tracing the geodesies, evaluating the local 
energy deposition rate [erg s""*^ cm~^] everywhere around 
the black hole and then integrating over volume outside the 
disc to obtain the net energy deposition rate Eiyi? (energy at 
infinity per unit time at infinity). 

The neutrino emission and annihilation is concentrated 
near the black hole, where the accretion time-scale is short 
and M may be assumed to be quasi-steady. Then Ei^i^ de- 
pends on four parameters that specify the steady disc model: 
accretion rate M, viscosity parameter a, mass of the black 
hole M, and spin of the black hole a. 
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2.1 Neutrino source: disc model 

As the disc matter spirals into the black hole, it is viscously 
heated: the gravitational energy is converted to heat. Part of 
the heat is lost to neutrino emission, and part is stored in the 
disc and distributed between nuclear matter, radiation, and 
pairs, in perfect thermodynamic and nuclear statistical 
equilibrium. The equilibrium microphysics is determined by 
only three parameters: temperature T, baryon mass density 
p, and electron fraction Ye (equal to the charged nucleon 
fraction). Other parameters — e.g. the electron chemical 
potential /ie and density of pairs n± — are derived from 
T, p and Ye. 

The neutrino emission peaks in the inner region of the 
disc. The far dominant emission mechanism is the e~/e^ 
capture onto protons/neutrons. 
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The neutrino-cooled discs are nearly in /^-equilibrium: the 
rates of the two reactions in equation ([T]) are practically 
equal, which determines the value of Ye{p^ T) (Imshennik et 
al. 1967; Beloborodov 2003). In principle, all three flavors of 
neutrinos could be emitted from the disc, but only electron 
neutrinos need to be considered; the emission rates for muon 
and tau neutrinos are negligible (CB07). 

A detailed model for neutrino-cooled relativistic discs 
was developed in CB07, and we use their model in our cal- 
culations. The approximate hydrodynamic disc equations 
are solved with the vertically-integrated a prescription. The 
equations include radial transport of heat and lepton num- 
ber. Local microphysics is treated exactly: nuclear compo- 
sition, electron degeneracy, neutrino emissivity and opac- 
ity etc., using the equilibrium distribution functions for all 
species except neutrinos. Neutrinos are modeled separately 
in the opaque and transparent zones of the disc, matching 
at the transition between the two zones. 

The disc model provides the vertically averaged T, p. 
Ye , pe (which are approximately equal to their values in the 
midplane — 7r/2) and the half-thickness of the disc H. 
The v and V spectra emitted from the neutrino-transparent 
(optically thin) zone of the disc are given by 
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for neutrinos (^ > 0), and 
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(3) 



for anti- neutrinos (E > ^ + 1). Here h = 27r/i, A = h/nieC, 
q = (rrin - mp)/me = 2.53, K = 6.5 x IQ-'^s'^ 6 = 
kT/nieC^, Pe is electron chemical potential in units of rrieC^ 
and E is neutrino/anti-neutrino energy in units of rrieC? . 

The spectra emerging from the opaque zone are con- 
trolled by neutrino transport through the disc, which cannot 
be reliably calculated — it depends on the unknown vertical 
distribution of viscous heating. Fortunately, the annihi- 
lation rate depends only on the energy fluxes Fi^ and F^? 



from the disc surface, which are insensitive to the neutrino- 
transport details. This fact has a simple analytical explana- 
tion (Beloborodov 2008). It follows from the proportionality 
aui) oc EuEu^ where a^u is the cross-section for annihilation 
for neutrino and anti-neutrinos with energies Ei, and E^^. 

To demonstrate that the rate of annihilation depends 
on Ft, and Ft, but not on the exact shapes of v and v spectra, 
we consider below two extreme models A and B, and they 
give practically the same annihilation rates: 

Model A: Neutrinos and anti-neutrinos are emitted 
with the same spectra as found inside the disc. The spectra 
are normalized so that the emerging emission carries away 
the known energy fluxes Fu and F^. In the opaque region, 
the spectra of v and v are described by Fermi-Dirac distri- 
butions. The temperature and chemical potential for v and 
V inside the disc are obtained from the numerical models of 
CB07. 

Model B: Neutrinos and anti- neutrinos are emitted 
with a thermal spectrum with zero chemical potential and 
temperature T = Teff , where Teff is the effective surface tem- 
perature defined by (7/8) crTgfj — F^ + F^ {a is the Stefan- 
Boltzmann constant and the coefficient 7/8 takes into ac- 
count the difference between statistics of photons and zy, u). 

When the disc is efficiently cooled (neutrino energy flux 
almost balances viscous heating), Teff is given by the stan- 
dard thin-disc model of Page & Thorne (1974): Teff ^ T|^. 
This regime occurs in a broad range of accretion rates 
Mign < M < Mtrap (CB07). If M < Mign, the disc tem- 
perature is not high enough to ignite the neutrino emit- 
ting reactions. If M > Mtrap, the emitted neutrinos become 
trapped in the disc and advected into the black hole. Our 
third (simplest) model for the neutrino source is defined as 
follows. 

Model C: Neutrinos and anti-neutrinos are emitted 
with a thermal spectrum that has zero chemical potential 
and the following temperature. 



Teff (M,r) = T|^ (Mign, r) 


X < (M/Mign)^/" 

(Mtrap/Mign)'/^ 



M < Mign 
Mign < M < Mtrap 
M > Mtrap 
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This model does not even require the calculation of the disc 
structure as T^^ is a known analytical function of r (Page & 
Thorne 1974). As we show below, this simplest model gives 
remarkably accurate result for Euu. 

The characteristic accretion rates Mign and Mtrap de- 
pend on the viscosity parameter a. They were calculated in 
CB07, and the numerical results are well approximated by 
the following formulae. 
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The coefficients i^ign and i^trap are functions of the black- 
hole spin a. For a = 0, K[gn = O.O71M0S~^ and iCtrap = 
9.3M0S"^ For a = 0.95, Ki^n = O.O21M0s"^ and Ktre^p = 
1.8M0S-^ 

In all three models A, B and C the neutrino emission 
is assumed to be isotropic in the local rest frame of the disc 
(which is in Keplerian rotation around the black hole). 
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2.2 Neutrino transport 



To evaluate the annihilation rate at a given point we need 
to know the local u and u distribution functions. They can 
be obtained using the known neutrino distribution functions 
at the surface of the disc. To a first approximation, neutrinos 
obey the collisionless Boltzmann equation, because most of 
them do not participate in any interactions (and eventually 
escape to infinity or get captured into the black hole). The 
Boltzmann equation in curved spacetime has the same form 
as in flat spacetime. It states that the phase-space density 
(or occupation number) of neutrinos remains constant along 
their trajectories. 



0. 



(6) 



Here x^(A) is a parameterized worldline for a or zy. The 
neutrinos emitted by the disc have huge energies compared 
to their rest mass and we treat them as massless particles 
propagating along null geodesies in Kerr spacetime. We use 
Boyer-Lindquist coordinates = (t, 0, r, 0), where the Kerr 
metric has the form. 



•dr^ + 
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The metric tensor Qa/s is specified by two parameters: rg = 
2GM/c^ and the spin parameter \a\ < 1; it is given in e.g. 
Chandrasekhar (1998). 

Equation ^ is covariant and its solution takes into 
account Doppler and gravitational redshifts. The null 
geodesies in Boyer-Lindquist coordinates are described by 
the known ordinary differential equations of first order (e.g. 
Chandrasekhar 1998) which we solve numerically. 

Figure [l] shows the image of the accretion disc observed 
from three locations near the black hole. Colour represents 
the redshift (or blueshift) of neutrinos as they propagate 
from the emission point on the disc to the observation point. 
The asymmetry of the images is caused by the rotation of 
the black hole. 



2.3 Local rates of creation 

2.3.1 Reaction uiy ^ e~ 

The rate of lyjy annihilation at a given point depends only 
on the local momentum distribution of u and We use 
the phase-space occupation number to describe this distri- 
bution; it is denoted by / for neutrinos and by f for anti- 
neutrinos. The local energy-momentum deposition rate is 
given by (see e.g. Birkl et. al. 2007), 



where 



A{p,p') 



-j^-j^Aip,p)[p +p 



f{p)f'{p), 



C P°P'°(1 
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where p and p' are the four-momenta of neutrino and anti- 
neutrino, is the angle between the spatial (3D) components 
of p and p, ao = AmlGl/irh^ ^ 1.71 x 10"^^cm^ Gf is 
Fermi constant, Ci = 0.78 and C2 = 1.06. We neglect the 
mass of ly and i.e. assume p^Pa — and p'^p^ — 0. 
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Figure 1. Projection of the sky for three different observers 
receiving neutrinos from a disc of radius 140r^. The colour codes 
the ratio of the energy at reception to the energy at emission for 
neutrinos, as measured in the local ZAMO frames (Section 2.4). 
All three observers are located at 3.1 horizon radii rh around a 
black hole of spin a = 0.95. Their positions are approximately 0, 
7r/4 and 7r/2. The observer sky is parameterized with two angles 
a e (0,7r) and /3 G (0, 27r) in the ZAMO frame. The axis a = 
points into the black hole and a = tt points away from the 
black hole. The x and y coordinates in the figure are {x, y) = 
(f sm/3,f cos/3). 
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2.3.2 Reaction v - 



in a strong magnetic field 



Code accuracy 



Neutrinos propagating in a strong magnetic field B can 
create electron-positron plasma through the process u 
z^e^e" (Kuznetsov & Mikheev 1997; Gvozdev & Ognev 
2001). A neutrino with four-momentum and energy E 
deposits energy-momentum into the plasma with the fol- 
lowing rate 
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where is the angle between the magnetic field and the 
neutrino momentum, Cv ~ 0.96 and Ca = 1/2. The local 
energy- momentum deposition rate via reactions ly z/e^e~ 
and u ue~^e~ is given by 



(11) 



2.4 Integration of the energy deposition rate over 
volume 



The local energy-momentum deposition rates (eqs.jsjand 11 ) 
can be calculated in any frame of reference. We chose the 
frame of Zero Angular Momentum Observer, ZAMO. These 
observers do not move in r or ^ directions. Their (j) mo- 
tion corresponds to zero angular momentum (p^ = 0). We 
use ZAMO frames because they are well defined everywhere 
outside the black-hole horizon (no observers can be static 
in Boyer-Lindquist coordinates inside the ergosphere). The 
local orthonormal coordinates of ZAMO, dx^ , are related to 
dx^ = (dt, dr, dO) by 



dx^ — A^i3 dx^ , 
where 





(12) 



det^ — \/—g and the element of four-dimensional volume 
measured by ZAMO dV is related to the Boyer-Lindquist 
coordinate volume dV — dt dcj) dr dO by 



dV ^ 
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Let dP^ — Q^dV be four- momentum deposited in dV as 
measured by ZAMO. The deposited energy measured by a 
distant observer is dE = —cdPt. It can be expressed in terms 
of dP^, 
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where we used the transformation dP^ — 
yields 
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The net energy deposition rate outside the horizon dE/dt 
(measured by a distant observer) is given by 



E = 



9tct> 



where ru is the horizon radius. 



-g d(\) dr dO^ 
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Half spatial resolution Twice geodesies 



a=0 



\/:^E/E\ = 0.03 \/\E/E\ = 5 X 10-^ 



a=0.95 \/:^E/E\ = 0.07 



|A£;/£;| = 0.02 



Table 1. Fractional change of the energy deposition rate when 
the grid resolution is reduced from 500 to 240 points (left column), 
and when the number of geodesies is increased from 5 x 10^ to 
10^ (right column). 



2.5 Numerical Method 

We compute the local rates of four-momentum deposition 

on a spatial grid. The problem is axially symmetric and 
the grid is set on the (r, 0) plane. It covers the region ru < 
T < Tmax and < ^ < 7r/2 (we use the symmetry about 
the equatorial plane — 7r/2). The radius rmax is chosen 
between 2Qrh and 38r/i, depending on the black hole spin. 
The grid has 25 x 20 points. Its spacing is logarithmic in the 
r-direction and uniform in the ^-direction. 

We calculate according to equations ([8| and (11). 
The distributions f{p) and f'{p) are Lorentz- invariant 
(scalar) functions of four-momentum. They remain constant 
along the neutrino trajectories (eq. |6|. Therefore, to calcu- 
late f{p) measured by ZAMO at a given point, it is sufficient 
to trace the neutrino with momentum p back to its emission 
point, find its four-momentum there in the rest-frame of the 
disc, pem, and use the equality f{p) = /em(pem)- The values 
of /em(Pem) and /em(Pem) are provldcd by the disc models 
described in § 2.1. 

For every point of the grid we trace back the neutrino 
trajectories for 5000 directions uniformly distributed on the 
local ZAMO sky. Each direction is followed until the geodesic 
reaches the disc, goes into the black-hole horizon or reaches a 
maximum radius that we set equal to 140r^. For trajectories 
coming from outside 140r^ or connecting to the horizon we 
set f — f' — ^- For trajectories connecting our grid point to 
the disc, we find / and f as described above. 

The accuracy of our calculation of can be estimated 
by doubling the number of sampled geodesies; it is a few per 
cent. The accuracy of the volume- integrated quantity E is 
controlled mainly by the number of grid points. We checked 
it by repeating the same calculation with a coarser grid; the 
estimated error in E is smaller than 10 per cent (Table 1). 



2.6 Comparison with previous works 

We compared our calculations with three previous works: 
Birkl et al. (2007), Asano & Fukuyama (2001) and Popham 
et al. (1999). 

Birkl et al. (2007) calculated Q^^ and the corresponding 
volume- integrated energy deposition rate E^u outside the er- 
gosphere for several toy models of the neutrino source. They 
traced exactly the neutrino trajectories in the Kerr metric. 
Unfortunately, they incorrectly computed Ejjf} by integrat- 
ing dP^ instead of —dPt. For test purposes we did a similar 
integration. We computed their models D and REF, which 
assume that neutrinos are emitted by an isothermal black- 
body ring (see Table 1 in Birkl et al. 2007). The results 
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Figure 2. Color-coded and contour plot for log^o Qtj^ around an 
accretion disc with M = IMq/s. The black hole has mass M = 
SMq and spin parameter a = 0.95. The horizon sphere (black) 
has the radius ~ l.SCM/c^, and the inner edge of the disc (the 
marginally stable orbit) is at rms ^ = 2GM/c^ . Arrows show 
the projection of Q]j^lQ\j^ on the plane of the figure. The white 
curve is where the radial component of injected momentum Q[j^ 
changes sign. It roughly indicates the region where the deposited 
energy may be lost into the black hole rather than escape in an 
outflow. Disc Model A with viscosity parameter o; = 0.1 (Section 
2.1) was used in the calculation. Practically the same ^ is found 
for Models B and C, and for different ol (e.g. ol — 0.01). 



agreed within the numerical errors of their and our calcula- 
tions. 

Asano & Fukuyama (2001) calculated E^j^j on the rota- 
tional axis in Kerr metric as a function of r. The neutrino 
source was modeled as a blackbody disc with the temper- 
ature having a power-law dependence on radius. For test 
purposes, we repeated the calculation for their isothermal 
disc model and obtained Euu{r) on the axis. The functional 
shape of Euu{r) agrees with the shape of function G{r) in 
Asano & Fukuyama (2001) (see their eq. 18). The normal- 
ization of G{r) appears to be arbitrary, so we were unable 
to compare the numerical values of Et^t^. 

Popham et al. (1999) made more realistic assumptions 
about the neutrino source, similar to our Model C (Sec- 
tion 2.1). They evaluated numerically Et^tj for several values 
of accretion rate M and black-hole spin a (see their Ta- 
ble 3). The neutrino geodesies in Kerr metric were replaced 
by straight lines. When comparing their results with ours 
(which are presented in the next section) we found a sig- 
nificant disagreement, exceeding factor of 10. Generally, we 
find that E^i? of Popham et al. (1999) was overestimated. 
The dependence of Ej^i, on a that is suggested by Table 3 in 
Popham et al. (1999) is incorrect, apparently steeper than 
what we find numerically and estimate analytically (see next 
section) . 
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Figure 3. Same as Fig.[2]but for a non-rotating black hole, a = 0. 



3 RESULTS 

The two processes of creation considered in this paper 
(Section 2.3) give the total energy deposition rate E = Euv + 
EuR. We first focus on hjui; and discuss the contribution E^,n 
separately in Section 3.2. 



3.1 Energy deposition from vv ^ e 

Figures 2 and 3 show Q^z^ for an accretion disc with M — 
IM0 s~^ around a black hole of mass M — 3Mq. In Figure 2 
the black hole is assumed to be rapidly rotating (a = 0.95), 
and in Figure 3 it is non-rotating (a = 0). In the case of a = 
0.95, the deposition rate Qti? peaks closer to the black hole 
and reaches much higher values, because the disc extends to 
smaller radii and emits a higher neutrino flux. 

Volume integration of Qtz^ and Qtu, as described by 



equation (16), gives the net energy deposition rate due to 
lyjy annihilation outside the black-hole horizon, Ejyj^. Figure^ 
shows Eiyu as a function of the disc accretion rate M for the 
two cases, a = and a — 0.95. For all M, Etjij is much higher 
when the black hole is rapidly rotating. 

The results are sensitive to M and a, but remarkably 
insensitive to the details of the disc model. 

(i) The uncertainty in the vertical structure of the ac- 
cretion disc leads to a small uncertainty in Euu as illustrated 
by two extreme models described in Section 2.1: Model A 
and Model B. The results of both models are well approx- 
imated by simple Model C (eq. [4|. Deviations of Model C 
from Model A arise mainly where Model C does not accu- 
rately predict the neutrino flux from the disc, i.e. where the 
disc is not strongly cooled by neutrino emission. 

(ii) The uncertainty in viscosity parameter a ^ 0.1 has 
almost no effect on Et^i, as long as Mign < M < Mtrap- In 
particular. Model C in this range of M is explicitly inde- 
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Figure 4. ^Jj^i^ as a function of M for a non-rotating black 
hole (a = 0) and for a rotating black hole (a = 0.95) of mass 
M = SMq . Open symbols show Model A and filled symbols show 
Model B (Section 2.1). The results of both models are well ap- 
proximated by simple Model C, which is shown in the figure by 
lines; the line is dotted at low M where the disc is transparent to 
neutrinos. 



pendent of a. The two characteristic accretion rates Mign 
and Mtrap depend on viscosity parameter a (see eq. [5|; in 
Figure [4] we assumed a = 0.1. 

The dependence of Euu on the black hole spin, for a 
fixed M = IM0/S, is shown in Figure [g] Instead of using 
the spin parameter a directly, it is more instructive to plot 
Euu versus radius of the last (marginally stable) orbit rms- 
Then one can see the power-law dependence of Eiy 1/ on Tms '• 
Ejyi} oc r~s"^. This power- law is accurate for rms > Tg which 
corresponds to a < 0.9. For rms < t^, Ej^i, has a somewhat 
stronger dependence on rms- The standard relation between 
rms and a (e.g. Page & Thorne 1974) is shown in Figure [s] 
For non-rotating black holes rms = QGM/c? and for maxi- 
mally rotating black holes rms = GM/c? . 



3.2 Energy deposition from v 
magnetic field 



ve^e in a strong 



The four- momentum deposition due to reaction v ve^e 



is given by equation (11 ). It depends on the magnetic field. 
We will assume the strongest field that could be expected, 

where Pmax is the maximum pressure in the disc; we find 
this pressure using the numerical disc model of CB07. For 
illustration, consider a toy model where the magnetic field is 
perpendicular to the disc, uniform in the region r < 25r^ and 
zero outside this region. The model overestimates a realistic 
magnetic field around an accretion disc. We will show that 
even such a strong field gives a modest Ej^b ■ The maximum 



Figure 5. Radius of the marginally stable orbit, rms, as a func- 
tion of the black-hole spin parameter a. 
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Figure 6. E^^i, and Ej^b as functions of rms for fixed accretion 
rate, M = 1Mq/s. 



pressure Pmax and B — (SyrPmax)^^^ depend on M and a. 
The magnetic field ranges in our models from ^GxlO^'^G 
to ^ 8 X 10^^ G, increasing with M and a. 

Figure [g] shows Ei^b as a function of rms for fixed M — 
1Mq/s. It approximately follows the same scaling relation 
as we found for ^zyzy(rms), Ei,b oc r~^"^. 

Figure [t] shows the ratio Eub/E^u as a function of M. 
This ratio is small for accretion rates above Mign and varies 
slowly with M or a. Thus, we find that E^^b makes a small 
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Figure 7. Ratio Ei^b /Eiyi? as a function of M for a = and a = 
0.95. The accretion disc is assumed to have viscosity parameter 
a = 0.1. 



Figure 8. Same as Fig. |4] but now including the contribution of 
reaction u — ^ h'e^e~. This contribution Ei^b is shown only at 
M < Mign5 where it is important; Ej^b is small compared with 
Ejyjy at M > Mign (cf. Fig, 
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contribution to the total E. This conclusion is valid for 
the most interesting range of accretion rates M > Mign- 
For M < Mign, the neutrino flux quickly decreases, which 
strongly suppresses the lyi? annihilation above the disc. The 
reduction in reaction u z/e^e~ is less severe and Eub ex- 
ceeds Ejyi?. As a result, Ei^b dominates the energy deposition 
rate E when M is below Mign. 

The dependence of the energy deposition rate on M 
is summarized in Figure [S] a modified version of Figure |4] 
It approximately represents the total E by showing Euu at 
M > Mign and Ej,b at M < Mign. 



3.3 Scaling of and Ej^b with M, M and rms 

The Ei^i^ dependence on M, M and a may be estimated 
using simple analytical arguments (Beloborodov 2008). The 
total rate of z/zy annihilation around the disc Ni^r, [s~^] is 
proportional to the volume of the main annihilation region 
r^, which scales as r^s? and to the typical local annihilation 
rate in this region, riuu acnuTiu. Here rii, and Uu are the 
number densities of neutrinos and anti-neutrinos and a is 
the annihilation cross-section. The cross-section scales with 
the energies of the annihilating v and v approximately as 
cr oc EjyEi^, which leads to 



(18) 



The energy fluxes Fi, and Fi) from an efficiently cooled disc 
(Mign < M < Mtrap) Carry the released gravitational energy 
and do not depend on details of the disc structure. The fluxes 
scale with M, M, and r as F^, oc MM/r^. This gives 



AT 3 . 

Nuu oc r rtuu oc — - 



(19) 



Here we assumed that the size of the main annihilation re- 
gion r is proportional to rms = XmsVg, where Xms is deter- 
mined by the black- hole spin parameter a (see Fig. 5). 

The net heating rate Ei^i? may be estimated as Ej^i) ^ 
{El, + Eu)Ni,u, where Ei^ ^ E^ oc F^^ is roughly estimated 
using Model C (Section 2.1). This gives 



Eiyiy OC T-^qFi 



9/4 



-15/4^9/4^ 



-3/2 



(20) 



Our numerical results confirm the scaling Et^i, oc M^^^ and 
give a somewhat steeper dependence on Xnis, tLih>h> CX 
(Section 3.1). The relation Xms(a) can be substituted here 
to get the dependence Ei^ip(M, M, a). 

It is instructive to compare Eub with Euu. From equa- 
tions ([8| and ([TT]) one can derive the following estimate. 



EuB ^ QIb 



'0.1af^\n(^E 



(21) 



where a/ = e'^ /he 1/137, Bq mlc^/eU ^ 4.4 x 10^^ G, 
E is the average energy of neutrinos in units of rrieC^, 
Uu Fu / CIS the neutrino energy density and Ub — i^^/Svr. 
One may expect that Ub is proportional to the pressure in 
the disc P, which depends on viscosity parameter a. Our 
numerical calculations gave Eus/Euv ^0.1 for discs with 
a = 0.1 and M > Mign (Fig. 7). Since P oc models 
with smaller a would give a higher E^^b oc However, 
for discs with M > Mign, it will not dominate over Et^i, in 
the plausible range of a > 0.01. 

For small accretion rates M < Mign , F^ is strongly sup- 
pressed and Euu is suppressed as F^^"^. Since Eub oc F^ its 
suppression is less severe. As a result Ei,b dominates over 
Ei^i^ when M < Mign as we found numerically in Section 3.2. 



8 Ivan Zalamea and Andrei M. Beloborodov 



4 CONCLUSIONS 

We have performed detailed numerical calculations of 
creation around hyper- accreting spinning black holes. We 
studied two reactions: e^e~ and f z/e^e~ (in 

a strong magnetic field). The reaction of annihilation 
dominates the energy deposition rate around discs with 
M > Mign, which are strong emitters of u and u. We found 
that the net energy deposition rate due to this process, Ei,i), 
is well approximated by a simple formula (see Fig. 4), 

C M<Mign ^ 

X } m^/* Mign<M<Mtrap >ergs-\(22) 
[ mj*p M>M,,,p J 

where m — M/Mq s""*", Xms = f^ms{o)jTg^ Tg — 2GM/c?^ 
and Mtrap, Mign are given in equation (Isl. The dependence 
of Et^z, on the black hole spin is huge: x~s" varies by a factor 
of 200 for < a < 0.95. Note that a (viscosity parameter 
of the disc) enters the result only through Mign and Mtrap- 
Our numerical simulations in this paper are limited to black 
holes with mass M = SM©, and the E^u dependence on M 
is evaluated analytically. 

The efficiency of annihilation can be defined as e = 
E^^jE where L is the total neutrino luminosity of the disc. 
For example, a — 0.95 (which corresponds to Xms = 0.97) 
gives L ^ 0.15Mc2 (CB07) and e ^ 0.05m^/^ for Mign < 
M < Mtrap. Note that Euu is defined in this paper as the 
total energy deposition rate outside the event horizon. A 
fraction of the created plasma falls into the black hole 
and does not contribute to the observed explosion (Fig. |2|. 
The corresponding refinement of e depends on the plasma 
dynamics outside the disc, which is affected by magnetic 
fields and hard to calculate without additional assumptions. 

The obtained Euu may be comparable to GRB luminosi- 
ties Lobs- A plausible typical value for Lobs is ^ 10^^ erg/s 
(it depends on the beaming angle of the observed explosion, 
which is usually hard to estimate from available data). This 
power is easily supplied by neutrino heating if the black hole 
has a large spin, e.g. a = 0.95, for a moderate accretion rate 
M > O.3M0 s"\ For a non-rotating black hole, this mecha- 
nism of GRB explosion requires 10 times higher accretion 
rates. 
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